Simple model squirmers with tunable velocity 
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We present a scheme of self-propelling liquid droplets which closely mimics the locomotion of some 
protozoal organisms, so-called squirmers. In contrast to other schemes proposed earlier, locomotion 
paths are not self-avoiding, since the effect of the squirmer on the surrounding medium is weak. 
Our results suggest that not only the velocity, but also the mode of operation (i.e., the spherical 
harmonics of the flow field) can be controlled by appropriate variation of parameters. 



PACS numbers: 47.63.mf, 47.20.Dr, 47.63.Gd 

Microscopic organisms have developed a variety of 
strategies to achieve locomotion in liquids at low 
Reynolds' numbers [T|. Flagellated cells such as certain 
bacteria, sperm, E.Coli, and trypanosomes use nonrecip- 
rocal motions of their flagella to overcome viscous forces 
and set themselves in motion [2]. Organisms such as 
cyanobacteria, Paramecium, and Volvox belong to a class 
of swimmers referred to as squirmers and are driven by 
tangential and/or radial deformations of the cell surface 
[3J. Squirming motion is very appealing for elucidating 
the hydrodynamics of microscale swimming, since the ve- 
locities in the near and far field around such a swimming 
organism can be well described analytically [U \5\ . This 
makes such swimmers ideally suited for the quantitative 
study of many open questions regarding the hydrody- 
namic effects on their interaction with surfaces and with 
each other, their behavior in external flow, and the ori- 
gins of coupled and collective behavior [3J. However, in 
contrast to the large number of natural micro-swimmers 
and analytical models, very few artificial self propelled 
objects have been realized [B], and even fewer mimic 
actual biological mechanisms |7_. Apart from potential 
technological benefits, a significant appeal of model sys- 
tems is the ability to provide quantitative control and 
minimize the individual variations pertinent to biologi- 
cal systems. 

We present here a simple model squirmer consisting of 
an aqueous droplet moving in an oil 'background' phase. 
Propulsion arises due to the spontaneous bromination of 
mono-olein (rac-Glycerol-l-Mono-oleate) as a surfactant, 
which is abundantly present in the oil phase, such that 
the droplet interface is covered by a dense surfactant 
monolayer. The bromine 'fuel' is supplied from inside 
the droplet, such that bromination proceeds mainly at 
the droplet surface. It results in saturation of the un- 
saturated C-C bond in the alkyl chain of the surfactant, 
thus rendering it a weaker surfactant. As we will see be- 
low, this results in a self-sustained bromination gradient 
along the drop surface, which propels the droplet due to 
Marangoni stresses. 

We demonstrate this squirmer scheme with nanoliter 
droplets containing 25 mM bromine water in a continu- 




Figure 1. Top: Path of a single squirmer droplet. The persis- 
tence length is clearly large compared to the droplet radius, 
indicating propelled motion (Scale bar: 300 microns). Bot- 
tom: Time lapse series of seven droplets in a mocrochannel. 
The droplets change direction without noticeable reduction in 
velocity. 



ous oil phase of squalane containing 50 mM mono-olein 
(MO). The critical micelle concentration (CMC) is 1.5 
mM. The droplets are confined by two hydrophobic glass 
plates to a quasi 2 dimensional space, thus simplifying 
droplet tracking. The top panel of Fig. [I] shows the tra- 
jectory of a single squirmer droplet over a duration of 
400 seconds. The velocity of the droplet in this duration 
is roughly constant at about 15 //m/sec. The trajec- 
tory is reminiscent of a random walk, with a persistence 
length which is larger than the droplet size and clearly 
far beyond what would be expected for Brownian motion. 
A particularly important observation is that the trajec- 
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tory crosses itself. This is in sharp contrast to other 
schemes, where the propulsion mechanism itself changes 
the surrounding medium strongly enough to prevent self- 
crossing of the path |5]. That this is clearly not the case 
here is demonstrated even more convincingly in the bot- 
tom panel, which shows a time lapse representation of 
seven droplets moving in a micro-channel. As two drops 
touch each other, they reverse their direction of motion 
and perambulate the channel again, without significant 
reduction in velocity. 
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Figure 2. Top: Schematic of a micro-droplet squirmer. 
Bromination increases the tension of the droplet surface from 
1.8 mN/m to 3.1 mN/m. The convective flow pattern (shown 
in the rest frame of the droplet) is accompanied by a gradient 
in the bromination density. The corresponding Marangoni 
stress propels the droplet. Bottom: velocity field around a 
droplet squirmer. The magnitude of the flow velocity (color 
code) and streamlines along a horizontal section through the 
center of a squirmer droplet are shown. Scale bar: 100 mi- 
crons; velocity scale (right) in microns per second. 

In order to gain some insight into the propulsion mech- 
anism, let us consider a spherical droplet with radius R. 
The total coverage, c, of the droplet surface with the 
mono-olein, either brominated or not. is assumed to be 



roughly constant and in equilibrium with the micellar 
phase in the oil. The brominated fractional coverage shall 
be called b. If the droplet moves, there is (in the rest 
frame of the droplet) an axisymmetric flow field, u{9), 
along its surface. The equation of motion for b is 



db 
di 



k(bo — b) + div (D^gradfo — ub) 



(1) 



where bo is the equilibrium coverage with brominated 
mono-olein (brMO). It is determined by the bromine sup- 
ply from inside the droplet and the rate constant, k, of 
escape of brMO into the oil phase. Di is the diffusivity 
of the surfactant within the interface. 

The droplet motion is accompanied by a flow pattern 
within the droplet and in the neighboring oil, which can 
be determined from u{9) The corresponding viscous 
tangential stress exerted on the drop surface must be bal- 
anced by the Marangoni stress, grad7(#) = Mgrad&(#), 
where 7 is the surface tension of the surfactant-laden 
oil/water interface, and M — dj/db is the Marangoni 
coefficient of the system. Expanding the bromination 
density in spherical harmonics, 



b{9) = V 6 m P m (cosf 



(2) 



we can express the velocity field |9] at the interface as 
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where C™ denote Gegenbauer polynomials, and /i is 
the sum of the liquid viscosities outside and inside the 
droplet. Inserting this into eq. and exploiting the or- 
thogonality relations of Gegenbauer and Legendre poly- 
nomials, we obtain 
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for all m > 0. We see that the different modes decouple, 
as far as linear stability is concerned. As long as b$M 
is small enough, the resting state is stable against fluc- 
tuations. However, when boM exceeds a critical value, 
the resting state is unstable, and the droplet sponta- 
neously starts to move. It is straightforward to see that 
for k < 3Di/R 2 , this happens first for the lowest mode at 
to = 1 which corresponds to a purely dipolar flow field. 

In order to determine the flow profile around a squirm- 
ing drop, we performed particle image velocimetry using 
a standard setup (ILA GmbH). The oil phase is seeded 
with 200 nm green fluorescent polystyrene beads (Duke 
Scientific) to be used as tracers. The result in the bottom 
panel of Fig. [2j shows that the flow profile indeed resem- 
bles a purely dipolar field. We note, however, that by 
careful adjustment of conditions, such as adding suitable 
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solvents to the oil phase, one should be able to tune k, 
and thus control which modes become unstable, thereby 
adjusting the flow pattern in the droplet and its vicinity 

To discuss the steady state velocity, V, we reconsider 
the total surfactant coverage, c. This adjusts itself as a 
balance between the molecular adsorption energy at the 
water/oil interface and the mutual repulsion of the ad- 
sorbed surfactant molecules. The equilibrium coverage 
thus represents a minimum in the interfacial energy. As 
a consequence, the interfacial tension will not change to 
first order if c is varied. Deviations of c from its equilib- 
rium value, which come about necessarily for any finite V, 
will thus be replenished from the surrounding oil phase 
without noticeable Marangoni stresses. 

For a dipolar flow pattern as shown in Fig. [2] we have 
divu(0) oc cos0. The density of surfactant thus takes the 
form p ps po + Spf(r) cos9, where Sp oc V. f(r) is some 
function of the radius. If D m is the diffusivity of the mi- 
celles, the diffusion current, D m dp/dr, must balance the 
depletion rate at the drop interface, cdivu(0) = 4^ cos# 
PJ]. As long as this holds, the only source of appreciable 
Marangoni stresses is the gradient in bromination den- 
sity, b(6). The drop thus keeps taking up speed accord- 
ing to eq. Q. This comes to an end when 8p ps p^. The 
surfactant layer at the leading end of the drop surface 
can then not be replenished anymore, and c comes sub- 
stantially below its equilibrium value. This leads to an 
increase of surface tension accompanied by a 'backward' 
Marangoni stress, and thus finally to a saturation of the 
velocity. According to the reasoning above, we expect 
that V « 2p D m /3c. 

There is no literature value for the diffusivity of MO 
micelles in squalane, but we can estimate it on the basis 
of the Stokes-Einstein relation assuming the radius of the 
micelles to be similar to the length of a MO molecule (2.3 
nm). Using 36 mPa.s for the viscosity of the squalane, we 
obtain D m = 2.6xl0~ 12 m 2 /s. We thus predict V/po ~ 

U - Z ' mM/1 ' 

This can be tested experimentally. In order to measure 
V, we used a reaction mixture within the droplets simi- 
lar to the Belousov-Zhabotinski (BZ) reaction, with reac- 
tant concentrations adjusted such as to prevent chemical 
oscillations. This results in a spatially and temporally 
constant bromine release rate in the aqueous phase for an 
extended period of time. Fig. [3^, shows the dependence of 
the droplet velocity as a function of the mono-olcin con- 
centration in the oil phase. The initial linear increase is 
in agreement with the above prediction (dashed line) . As 
the surfactant density is further increased (and thus the 
velocity), the complex exchange processes between the 
water/oil interface and the micelles will finally become 
the rate limiting step. As a consequence, the velocity is 
expected to level off, in agreement with our data. 

Our model predicts the locomotion velocity to be inde- 
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Figure 3. Squirmer velocity in units of microns per second, as 
a function of (a) surfactant concentration, where each point is 
an average of 50 different squirmer droplets each of diameter 
~ 80 microns and (b) time when an oscillating chemical reac- 
tion (BZ) takes place in the droplet. For the latter case, the 
optical transmission of the droplet is plotted in red (arbitrary 
units, linear scale) along with the velocity trace. 



pendent of the bromine release rate. This can of course 
not be strictly true in general, and it is instructive to 
demonstrate this experimentally. We take advantage 
again of the BZ reaction, this time creating an oscil- 
lating bromine concentration inside the droplet. This 
oscillation can be easily visualized by the optical trans- 
mission of the droplet. Care was taken that there were 
no spatio-temporal patterns within the droplets, in con- 
trast to some other recent work on droplet locomotion 
[TO] . Fig. [3Jd shows an overlay of the the BZ chemical os- 
cillations and the velocity of the squirmer droplet, which 
are anti-phase with each other. As it is known for the 
BZ reaction, the decrease of the transmitted intensity 
corresponds to an increase in the bromine concentration 
within the droplet. It is obviously possible to control 
the droplet velocity in some range, and we have demon- 
strated here how this can be done even in an autonomous 
manner. We expect this model squirmer scheme to be 
instrumental in elucidating the hydrodynamics of collec- 
tive behaviour, in particular concerning locomoting or- 
ganisms. 
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